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ABSTRACT 

We compared optical spectroscopic and photometric data for 18 AGN galaxies over 2 to 
3 epochs, with time intervals of typically 5 to 10 years. We used the Multi-Object Dou¬ 
ble Spectrograph (MODS) at the Large Binocular Telescope (LBT) and compared the 
spectra to data taken from the SDSS database and the literature. We find variations in the 
forbidden oxygen lines as well as in the hydrogen recombination lines of these sources. 
For 4 of the sources we find that, within the calibration uncertainties, the variations in 
continuum and line spectra of the sources are very small. We argue that it is mainly the 
difference in black hole mass between the samples that is responsible for the different de¬ 
gree of continuum variability. In addition we find that for an otherwise constant accretion 
rate the total line variability (dominated by the narrow line contributions) reverberates 

3 

the continuum variability with a dependency AL;,,,,, oc (ALcont.)-- Since this dependency 
is prominently expressed in the narrow line emission it implies that the luminosity dom¬ 
inating part of the narrow line region must be very compact with a size of the order of 
at least 10 light years. A comparison to literature data shows that these findings describe 
the variability characteristics of a total of 61 broad and narrow line sources. 

Key words: Galaxies: Active - Variability. 


1 INTRODUCTION 

Seyfert galaxies and quasars are the main types of active 
galactic nuclei (AGNs) in which the nuclei are dominated by a 
very powerful electromagnetic emitter. Seyfert nuclei are 10 to 
100 times brighter than the host galaxies they reside in. Thus, 
these nuclei are responsible for most activity that occurs within 
the central 10-100 pc of these hosts. 

Seyfert galaxies can show a wide variety of forbidden 
and permitted emission lines, which can have widths up to ~ 
10 OOOkms^* via Doppler broadening. Depending on the width 
of the emission lines, one distinguishes between two main cate¬ 
gories, Seyfert-1 (SI) and Seyfert-2 (S2): S2 galaxies show only 
narrow lines (< 1000kms“'), while SI show additional broad 
components (£ 10 OOOkms^*) in the permitted lines. The emis¬ 
sion lines are commonly believed to stem from the “broad line 
region” (BLR) and the “narrow line region” (NLR), which are 
gas clouds orbiting a central super-massive black hole (SMBH) 
at very high velocities. The gas clouds are excited by contin¬ 
uum radiation produced by an accretion disk of material infalling 
around the central SMBH. According to t he unified model (e.g., 
lAntonuccH 1 19931 : llJrrv & Padovanil Il99.5h . the existence of SI 
galaxies which have broad and narrow lines, and S2 galaxies 
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which only have narrow lines, can be attributed to inclination ef¬ 
fects. One possible scenario is a dusty torus which surrounds the 
accretion disk and BLR, obscuring the emission from the BLR if 
the galaxy is observed with a viewing angle parallel to the torus 
plane. 

Narrow-line Seyfert-1 (NLSl) galaxies are SI galaxies 
that have broad components with comparably narrow widths 
(< 2000 km s“*). From an observational point of view it appears 
that there are two main differences between S2 and NLSl 
galaxies. The first one is that NLSl objects have Fe ii in their 
spectra ( which originate in the BLR ), which is absent in the S2 
spectra IGiannuzzo & StirDall996h . The second is that some 
NLS1 sources show strong lines of highly ionized i ron in their 
spect r a e.g. [Fe viil/i5721 A and [Fe x]/l6375 A dOsterbrockl 
Il985l : lOsterbrock & Poggeiri98.‘5h . these lines are rare and not 
typical for the S2 spectra. 

AG Ns are characte rized by variability at almost all wave¬ 
lengths (IPetersonll200lh . Studying the variability on all time 
scales is a very useful approach to investigate and understand 
the physical properties of AGNs, there are many systematic 
studies on the variability of AGN using spectroscopic and 
photometric observations. Photometric measurements have 
been used to study the variability, despite the fact that bright 
line emission is located within the corresponding passband 
dVanden Berk et al.l2004h . Studying the spectral variability both 
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in line and in continuum emission has become an important 
topic that may help to outline differences between different 
Seyfert source classes. Variability may indicate how the active 
galactic nuclei (AGN) affect the regions surrounding them 
(i.e. the BLR and NLR) and how these regions change as a 
function of time. Many interesting results have been obtained 
from AGN samples of various types and sizes. A noted one is 
the anti-correlation between the conti nuum luminosity and vari - 
ability, which was first reported by lAneione & SmithI (Il972h . 
Later addit i onal i nformation on th i s effe ct was provided by 


iHooketaf Jl994l); Zuo et al d2012 

. 20131). Additionally, other 

authors dCid Fernandes & Sodrel 

200Cl: iMeusinger & WeissI 


l2013h have clarified that this anti-correlation can qualitatively 
be described via the standard accretion disk model, proposing 
that the variability is caused via the variation of the accretion 
rate. 


Since spectroscopic studies can be very time consuming, 
there are only a few detailed, spectroscopy-based comparative 
studies on the variability of AGNs. They either concentrate 
on reverbera tion mapping an alysis of extensively observed 
sources te.g.. lKasDi et al.ir2007h or d epend on just a few epochs 
of spectroscopic observations (e.g.. IWilhite et al]l2005h . Multi 
epoch spectroscopy observations have advantages compared 
to photometric observations, as they allow us to study the 
continuum and line emission at the same time. Excepting 
regions of high emission line density, these observations also 
allow us to extract the underlying continuum spectrum over a 
broad wavelength range and to study the spectral shape of the 
lines and the line spectrum in general. The emission line spectra 
then allow us to study the BLR and NLR regions surrounding 
the compact variable continuum nuclei. 


Examples of extensive investigations on the variability in 
Seyferts and quasars concentrating on the continuum and th e 
fluxes of the emis si on lin e s are, e.g., PetersoneLaL ( 198^ 
IPeterson & Gaskel] ( Il986h : IPeterson et al. (Il998h : IGuo & Gul 
(l2014h . While the BLR lines are found to reverberate the con¬ 
tinuum variability on time scales of days, there are also a few 
reports of long term line va riability in the NLR region: B ased on 
a time coverage of 8 years Iciavel & Wamstekeil (Il987h find for 
3C390.3 that the NLR is photoionized and reverberates the con¬ 
tinuum flux. They set an upper limit of 10 light ye ars for the size 
of the NLR. Based on a time coverage of 5 years IPeterson et ^ 
( l2013h find for the SI galaxy NGC 5548 that the [O III] AA4959, 
5007 emission-line flux varies with time and give a size esti¬ 
mate of 1-3 pc (i.e. 3-9 light years). There is a clear connec- 
tion between the con tinuum and line variability. Eor instance, 
IPeterson et al.l j 19841) found that for ~73% of the cases they 
studied, there is variability in the H/3 line flux corresponding to 
the changes in the continuum. On the other hand, in the same 
project, they found that this correlation is not true for the sources 
if the continuum changes by more than ~ 70%. They suggest 
that for cases in which the Hj3 flux does not correspond to the 
continuum, light travel time effects may have to be included. 
Furthermore, there is still no absolute conclusion on the corre¬ 
lation between the emission lines width and the line luminos- 
ity. An anti-correlation was found in a sample of 85 quasars by 
iBrotherton et all jl994|), with single-epo ch spectroscopic data. 
On the other hand Wilhite et~^ ( l2005l) found that there is a 
correlation between the emission lines width and the line lumi¬ 
nosity. They conclude this from a spectroscopic sample of 315 
quasars. 


2 PAPER OUTLINE AND SAMPLE SELECTION 

In this paper we investigate the continuum and narrow line 
variability of a small sample of objects over a good fraction of a 
decade. The spectra cover a major fraction of the entire optical 
wavelength domain from the blue to the red, including several 
prominent lines, and also have an appreciable signal to noise ra¬ 
tio in the continuum emission. For such a study it is essential 
that the spectra are well calibrated such that they can be com¬ 
pared to each other. It is also important to correct for iron emis¬ 
sion in order to extract proper values for the line luminosities of 
species other than iron. As the data has been taken with different 
instruments and different effective apertures on the sky it is fur¬ 
thermore essential to make sure that aperture effects, which may 
occur due to slit losses, or from extended and spatially resolved 
line or continuum emission from the host galaxies, do not con¬ 
taminate the results of the variability assessment. In the literature 
one rarely finds appropriately calibrated or sufficiently described 
data sets that are suitable for long term variability studies. Hence 
we made the effort to collect and combine a representative data 
set. This comprises a total of 18 sources (see summary in Sec. 
16. H . For 8 objects we obtained and/or performed the detailed 
analysis over the available optical wavelength range using data 
from the LBT and other observatories (see Sect. l3.1l and Tab.[T). 
These 8 sources were selected on the basis of their observabil¬ 
ity (given the allocated observing dates), their brightness (such 
that with the chosen integration times a decent signal to noise 
could be reached that allows comparison to public survey data), 
and the fact that a previous first epoch spectrum had been taken 
typically 5-10 years ago. Based on their [Mgll] absorption, the 8 
sources include 3 loBAL QSOs. However, given that their spec¬ 
tra also contain prominent broad emission lines, we decided to 
include them in the characterization of the BLQSOs/Syl vari¬ 
ability. In order to reach our goal of studying the difference in 
the variability properties of BLQSOs/Syl and NLSyl/S2 galax¬ 
ies, it was clear that literature data on additional sources had to 
be added to enlarge the statistical basis. Hence, all conclusions 
drawn in this paper are based on the resulting larger samples of 
sources. 

For 5 sources (3 NLSl and only 2 BLSl) we obtained spec¬ 
troscopic information from the literature. For a further 5 sources 
(2 NLSl and 3 BLSl) we found sparser but suitable line and 
continuum information in the literature. A detailed description 
of the observations and data reduction is given in Sect. fol¬ 
lowed by a description of the methods used for our analysis in 
Sect.|4] In Sect. |5] we present the results for all galaxies that we 
used in our study. 

The sample also comprises sources with a range of iden¬ 
tifications. For the purpose of this investigation we combined 
sources that show very prominent broad emission lines, i.e. 
Broad Line Seyfert 1 and Quasi Stellar Objects (in the following 
BLSl and QSO) in addition to the narrow lines on which our 
investigation focuses. We also put Seyfert 2 and Narrow Line 
Seyfert 1 (sf]and NLSl) in one group as their spectra are dom¬ 
inated by narrow line emission. An additional justification for 
this combination is based on the fact that while the unified model 
of Seyfert galaxies suggests that there are hidden broad-line re¬ 
gions (HBLRs) in all Seyfert 2 galaxies, there is increasing evi¬ 
dence for the presence of a subclass of Seyfert 2 sourc es lacking 
these HBLRs dZhang & WandEoOfil : iHaas et alj|2007l) . In these 
cases one cannot exclude that we have a free sight to the nucleus 
despite of the fact that the source is classified as a S2 galaxy. We 
also refer here to the discussion in the review on NLS1 galaxies 


* One of our sample sources is a composite HII/S2 galaxy - but see 
discussion in Sect. E] 
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by iKomossj ( 1200 Sh . Hence, if one classifies sources based on 
the presence of absence of pronounced broad line emission, then 
S2 and NLSl may be more comparable to each other than e.g. 
BLS1 with NLS1 sources. In the discussion in Sect.l6land l6.ll we 
analyze the observed degrees of continuum and narrow line vari¬ 
ability. In retrospect we found that the above described combi¬ 
nation of sources is indeed justified by their dilferent variability 
characteristics. 

Clearly, at some point such an investigation needs to be per¬ 
formed using larger samples. Since an appropriate baseline for 
variability studies is several years, such an effort needs time and 
can only be provided in the near future. However, to study first 
order elFects, the number of sources needs to be only sufficiently 
large to separate the median or mean properties within the sta¬ 
tistical uncertainties. As we outline in the paper, this can already 
be done with the current, representative sample presented here. 
In fact, the interpretation we present in Sect. 16.21 may also be 
taken as a prediction of the effect that the degree of the contin¬ 
uum variability is indeed reverberated by the degree of narrow 
line variability. This prediction was motivated by the results of 
the analysis of the small but representative sample we present in 
this paper. A conclusion is given in Sect.|7] 


resolution R of the instrument is about R ~ 2000, covering a 
wavelength range between 3800 and 9200 A. 

In this research we used the SDSS active galaxies (Seyfert 
&. Qu asar) Catalogs based o n the seventh data release of the 
SDSS dAbazaiian et al.ll2009h . 


All our project targets (except J035409.48-t024930.7 
and J015328-I-260939 which are not covered by the SDSS 
survey) can be found in different SD SS catalogs re- 


leases (DR5, DR6, DR7, DR8. and DR9) dYork et al. 

200C; 

Stoushton et alJl2002l: Cibson et al 

20091: Skrutskie et al. 

2006; 

Adelman-McCarthv et alj 

20081: 

Adelman-McCarthv & et al. 

2009: Abazaiian et alJl200^ 

: Adelman-McCarthv & et al. 2011: 

Inada et al.ll2012h. 


Typical magnitudes of SDSS sources e.g. in the g-band are 
in the range of 14.47 < g < 19.21, the FWHM of the combined 
bright narrow and broad lines typically range from about 300 to 
3000 km s“'. 

In this paper we analyze the optical spectra of eight galax¬ 
ies, one is classified as composite HII/S2, 3 sources are classified 
as SI (BLSl & NLSl), and 4 sources are classified as QSOs. 
Further information about the analyzed objects is presented in 
Tab.[T] In Figs.l^and lAll we show images of the sources and a 
nearby reference star extracted from the same image frame. 


3 OBSERVATION AND DATA REDUCTION 

In this section, we describe the telescopes and instruments 
we used for the spectroscopic observations and the procedures 
followed for the data reduction. We also use complementary data 
from the literature and public archives li ke the seventh data re¬ 
lease of the SDSS dAbazaiian et al.ll2009l) . 

3.1 Multi-Object Double Spectrograph (MODS) 

For our observations we used th e MODS spectrograph a t 
the Large Binocular Telescope (LBT) dPogge et al.ll20ld . l2012h . 
The LBT is located on Mount Graham in the Pinaleno Moun¬ 
tains southeastern Arizona, USA. The MODS spectrograph con¬ 
sists of a pair of identical double beam blue & red optimized 
optical spectrographs (MODS 1 &. MODS2) which work as indi¬ 
vidual spectrographs. The instrument can be used for imaging, 
long-slit, and multi-object spectroscopy. The charge-coupled de¬ 
vice (CCD) detector of MODS is an e2vCCD231-68 Skx3k with 
a 15 yum pixel pitch. 

We use MODS in long-slit mode with a slit width of B.'S. 
The spectrograph operates between 3200-10500 A with a spec¬ 
tral resolution of R=/l/(5/l ~ 2000 (i.e., ~ 150 km s“*) where the 
full spectral range of the grating spectroscopy is split into a blue 
and red channel, which are 3200-6000 A (blue) and 5000-10500 
A (red). All objects which we targeted with our LBT observing 
program are listed in Tab.[T] 

3.2 The Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (SDSS) is a sensitive pho¬ 
tometric and spectroscopic public survey which started sci¬ 
ence operation in May 2000 and cover s 10** square degrees o f 
the celestial sphere in the northern sky dStoughton et al.ll2002h . 
The survey uses a dedicated 2.5-m wide-angle optical telescope 
mounted at Apache Point Observatory in New Mexico, United 
States (Latitude 32‘’46'49'.'30 N, Longitude 105°49' 13750 W, El¬ 
evation 2788m). This instrument provides images and photomet¬ 
ric parameters in five bands (u, g, r, i, and z) with an average 
seeing of T.'5 and down to a limiting magnitude of ~ 22.2 in r 
band. Spectroscopy of selected objects is done through 3"fibers 
with calibration uncertainties of the order of 2%. The spectral 


3.3 Data reduction, wavelength calibration, and flux 
calibration 

Here we describe the data reduction of our spectroscopic 
observations, wavelength calibration, flux calibration, and the 
complementary data from public archives that has been used in 
our analysis. 

For reducing the data of the sources listed in Tab.[T] we used 
the python software modsCCDRed, which is provided by the 
instrument team. This software package gives us the possibility 
to create bias correction, flat fielding, normalization, and other 
standard steps of data reduction. 

First we create normalized spectral flat field frames for each 
channel (blue and red) through the following steps: 

• bias correction of the flat fields images; 

• median combination of the bias-corrected flats; 

• interpolation of the bad columns using the bad pixel lists 
for the detector; 

• elimination of the (wavelength dependant) color term to 
produce a normalized “pixel flat”. 

More information on producing the normalized calibra¬ 
tion frames can be obtained from the instrument manually Af¬ 
ter generating the flat, all bias-corrected science data were flat 
fielded. The science frames were then wavelength-calibrated 
using a wavelength-calibration map generated from lamp (ar¬ 
gon, neon, xenon and krypton) data, as provided by the in¬ 
strument. In this step, we used the reduction package Iraf 
to identify the lamp lines and then transformed the science 
frames accordingly. Furthermore, we tested the wavelength- 
calibration by using an OH-sky line at las dOsterbrock et al.l 
1 19961 ; lOsterbrock. Fulbrieht & Bidil 19971) and found that both 
calibration methods (skylines and lamp) are in good agreement 
with an uncertainty < 3A. 

For the next step, we flux-calibrated the one-dimensional 
spectra making use of two standard stars: G191B2B and Feige67 
for which data was taken under 1.3” seeing conditions. Both 
standard stars result in a consistent and independent calibration. 

For the final data product we used both stars for calibration. We 

^ http: //WWW . astronomy. Ohio- state. edu/MODS/Manuals/MODSCCDRed. pdf 
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wavelength [A] 


a 



wavelength [A] 


Figure 1. The standard stars Feige67 and G191B2B. 


used reduction routines on two different software platforms: IDL 
and IRAK We extract the one-dimensional spectra from the two- 
dimensional science frames of the two calibration star using rou¬ 
tines in the Iraf package. For test purposes we calibrated the 
two stars with each other and compared the results with pub¬ 
lished data. We also downloaded spectra of these two standard 
stars from the European Southern Observatory 0 and compared 
them with our results (see Fig.[T). Calibrating the science objects 
with the two reference stars independently resulted in spectra 
that agreed very well with each other within an average uncer¬ 
tainty of 2% and 3% in the blue and red channel, respectively. 
We performed that calibration for each source for the red and 
blue channel and combine these spectra to a single spectrum per 
source using Iraf. In this process, we found that there are no 
significant continuum offsets in the overlap area between the red 
and blue side. While this conhrms the consistency of our cal¬ 
ibration, the uncertainties in comparing data between different 
telescope s and instruments is more of the order a few percent, 
(see e.g.. lPeterson et al1l2002h i.e. of the same order as reached 
for the SDSS data releases. 

In order to ensure a sufficient inter-calibration accuracy be¬ 
tween telescopes and instruments we applied a correction for 


^ http: //WWW. eso.org/sci/observing/tools/standards/spectra/ 


slit losses. We observed our target sources (including the two 
calibrated stars; see above) in typically 1.2 to 1.4 arcsec seeing 
through a 0.8 arcsec slit. For the two extreme cases we applied 
a seeing correction based on Fig. lA2l For J0153 taken in 1.90” 
seeing we scaled the LET fluxes up by a factor of 1.35. For J1203 
taken in 0.82” seeing we scaled the LET fluxes down by a factor 
of 0.70 SDSS. 

In order to correct for a time variable slit loss contribution 
we proceeded in the following way: For each target we have 2 to 
5 exposures such that we have a good statistical estimate on the 
flux loss due to the combination of variable seeing and misplace¬ 
ment of the slit. In Fig. lA2l we show a number histogram of the 
intensity drop measured in the [OIII] 45007, [Oil] 43727, Y{f} or 
Hff line in the corresponding red and blue channel for the fainter 
exposures with respect to the brightest once. To first order we 
correct for this effect by adjusting all exposures to the flux level 
of the brightest exposure per source. Under the assumption that 
for the brightest exposure the slit was always centered on the 
source no additional correction for slit losses is necessary (i.e. a 
factor of unity; case a in Fig. lA2t . However, if we assume that 
no significant systematical misplacements happened and that the 
statistics for slit losses for the brightest exposures is similar to 
the statistics of the fainter exposures then an unfavorable and un¬ 
likely case will be that all brightest exposures were taken under 
conditions of a mean slit loss (case p in Fig. lA2t . For a final sec¬ 
ond order correction for slit losses we assumed a case between 
case a and case p and chose the mean between unity and the 
factor for mean losses of 1.17±0.10, i.e. a second order correc¬ 
tion factor of 1.08, after having applied the first order correction. 
Judging from Fig. lA2l the uncertainty on this correction for the 
time variable slit losses is probably of the order of 5% for the 
1.2” to 1.4” seeing cases. For J0153 with 1.9” this correction 
is probably overestimated by <5% For J1203 with 0.82” seeing 
our correction may be underestimated by 5-10% especially if 
one takes into account that for a better seeing slit positioning can 
also be done more reliably. Given the seeing and aperture combi¬ 
nations for the Eeijing and Hiltner telescope measurements (see 
section lSTl and lS]^ as well as the high SDSS calibration quality 
(see section no seeing/slit loss corrections were applied to 
these data. 

As we flux calibrate our observed spectra with the flux cal¬ 
ibration stars following the procedure presented above and show 
that the galaxies are all very compact (see section 3.2 and Fig.l 
and A.2 and and section ISTI for the marginally extended source 
J0938) we assume the inter-calibration uncertainties to be less 
than 10%. 


3.4 Magnitude measurements 

In the SDSS survey magnitudes are derived using different 
methods (Model, Cmodel, Petrosian, and PSF). For our work we 
used magnitudes derived via the PSF method in which a Gaus¬ 
sian model of the PSF is fitted to the object (dominated by a 
compact PSF like source). The method also accounts for the 
variation of the PSF across the field and was hence considered 
as most suitable for our study. 

We compared the flux density value of the photometric ob¬ 
servation^ with those obtained from spectroscopic observation 
for each source. In this way we can get an estimate on the degree 
of variability that is present in the nuclei of these objects. 

In our LET observations we used a long slit with aperture of 


^ We use the Gemini flux density/magnitude converter 
(http://www.gemini . edu/?q=node/11119l to convert the magni¬ 
tudes derived from the photometric images in the u,g,r,i, and z-bands 
into physical units (flux density F,i in erg s“* cm“^ A“*). 
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O'.'8, while for the SDSS observations a spectroscopic hber with 
radius 3" was used. To demonstrate that aperture effects between 
slit and hber measurements are negligible, we show SDSS z- 
band images of the galaxies together with images of nearby stars 
from the same frame (Figs. l^and lATt . We scale the stars to the 
Oux level of the galaxy and subtract the frames from each other. 
In all eight cases, no emission is left. This shows that the sources 
are all dominated by emission from an unresolved combination 
of a stellar bulge and an AGN. 
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Table 1. Coordinates, redshifts, classification and observing parameters for all sources. 


o^ 


name 

RA 

Dec. 

Redshift (z) 

Classi. 

Obs. 

date 1 

Obs. date 2 

Photo. 

obs. date 

Dl 

Ref. 


(hh mm ss) 

(O , 



Telesc. 

Da.(M Y) 

Da.(M Y) 

Da.(M Y) 

Mpc 


SDSS J120300.19+162443.8 

12 03 00.1 

16 24 43.8 

0.16552 + 0.00001 

NLSl 

SDSS 

04 2007 

02 2012 

06 2005 

794.5 

1 

SDSS J093801.63+135317.0 

09 38 01.6 

13 53 17.0 

0.10056 + 0.00001 

HII/S2 

SDSS 

12 2006 

02 2012 

01 2006 

463.2 

2 

SDSS J034740.18+010514.0 

03 47 40.1 

01 05 14.0 

0.03149 + 0.00003 

NLSl 

OHP" 

09 2003 

02 2012 

11 2001 

138.1 

3,4 

SDSS J115816.72+132624.1 

11 58 16.7 

13 26 24.1 

0.43966 + 0.00035 

QSO 

SDSS 

05 2004 

01 2012 

03 2003 

2429.5 

5,6 

SDSS J080248.18+551328.9 

08 02 48.1 

55 13 28.8 

0.66406 + 0.00035 

QSO 

SDSS 

01 2005 

01 2012 

11 2003 

3993.9 

5,7 

SDSS J091146.06+403501.0 

09 11 46.0 

40 35 01.0 

0.44121 + 0.00099 

QSO 

SDSS 

01 2003 

02 2012 

12 2001 

2439.7 

5,6 

2MASX J03540948+0249307* 

03 54 09.4 

02 49 30.7 

0.03600 + 0.00008 

BLSl 

Hiltner* 

10 2004 

02 2012 

— 

158.4 

3,8 

GALEXASC J015328.23+260938.5* 

01 53 28.2 

26 09 39.1 

0.32640 + 0.00127 

QSO 

Beijing‘S 

02 2002 

01 2012 

— 

1711.9 

9,10 


Refere nces for redshifts and po sitio ns: (I )^delman-McCarthy & et al (2)\Sdnche^^ Almeida et al\ l20I^: (3 )\Veron- Cet^ Vero}i\200di) : (4) \Hewitt & Burbide^ il7997l) ; (5} \Hewett & Wild l2Q7d) ; (6) \Zhans et al\ 

tlOldl) : (7) \Gibson et alWlOO^} : (8) \Chanan^owne^^^arffoiiU9^ ): {9) \Darlins ii^Giova^lli \200A : (}0) \zhou et al.\i200i) . 

Notes: The 3 QSOs (J1158, J0911, J0802) contained in the SDSS catalogue have been classified as low-ionization broad absorption-line (loBAL) due to a prominent [Mg ii] absorption in their spectra. The second spectroscopy 
obserx’ation have always taken by LBT The photometric obserx’ations were always taken by SDSS. For the final two sources in the table no SDSS photometric obserx’ations are available. 

* Not covered by SDSS surx’ey 

" Telescope of Obserx’atoire de Haute-Provence (OHP) 

^ Hiltner Telescope 
^ Beijing Obserx’atory 


Y. E. Rashed et al. 























Variability in Active Galaxies 1 


Furthermore, we fitted Gaussian functions to the images of 
the galaxies and the stars. The results are listed in Tab. [4] sup¬ 
porting the finding that the galaxy emission is dominated by the 
unresolved point source. In Tab. [2] we list the seeing values for 
the different exposures. 


4 ANALYSIS 

We developed a Python routine to manually fit the stellar 
continuum, the power-law contribution from the AGN, as well as 
the Fe ii emission and subtract all of them. In the residual spec¬ 
trum we can then fit the non-Fe ii emission lines. In following 
we explain each step: 


4.1 Stellar continuum subtraction 



restframe wavelength [A] 



restframe wavelength [A] 


To obtain an accurate measurement of the nuclear emis¬ 
sion line fluxes and the equivalent widths, we have to subtract 
the stellar component of the host or its bulge component. To re¬ 
move the stellar component we fitted a stellar population synthe¬ 
sis model to the entire spectrum. The templates used are given 
by a sample spectrum built by a population synthesis routine 
in iBruzual & Chariot ll2003h . The template for the young stel¬ 
lar population spectrum is taken 290 Myr after a starburst with 
steady star formation rate over 0.1 Gyr and solar metalicity. For 
the intermediate-age population we used template of a 1.4 Gyr 
old simple stellar population with solar metalicity. 

In Tab.|5]we show the flux contributions of the stellar com¬ 
ponent, the AGN/powerlaw component, and the Fe ii template, 
as well as the extinction Ay of the stellar component. We mea¬ 
sured the fraction of each component compared to the total flux 
in the (restframe) wavelength interval 5100 A to 5600 A. The 
main purpose of the subtraction was to remove the Fe ii emis¬ 
sion lines from the spectrum, particularly those blended with 
other emission lines. Since in most galaxies there are no obvi¬ 
ous stellar features, it is difficult to distinguish between different 
stellar populations. Therefore, we sum up the contributions of 
the intermediate-age and the young stellar population. Also the 
fit of the extinction Aj/ is not reliable, but was only chosen to 
fit the continuum slope as accurate as possible. The data for the 
line variability as discussed in Sect. includes all corrections 
derived from the fitting described here. 


4.2 Fe II subtraction 


The optical emission-lines of Fe ii are typical features of 
Seyfert 1 galaxies and quasars. There is huge variety in the 
strength of these lines, from very clearly shaped and detectable 
lines in some AGN spectra, all the way to very weak and imper¬ 
ceptible lines in other sources. iBoroson & GreenI ( Il992h (in the 
following BG92) reported a strong inverse correlation between 
the strength of Fe ii emission-lines and the width of the broad 
H/I line as well as the strength o f the [O III] line. This result was 
confirmed by other groups (e.g., Cqrbiii|[l997 ; Zhou et aT]| 2006 l : 


iKovacevic. Popovic ife DimitriieG d2010l : lDong et alJ201lh . For 

our work we used the Fe ii template from BG92, fitted and scaled 
to our galaxy spectra as shown in Fig. [3] Tab. |3 shows that the 
contribution of Fe ii in the given wavelength range is up to 23% 
and a proper subtraction of Fe ii is required, so the fitted spec¬ 
tra were then subtracted from our spectra. Finally, we obtained 
a spectrum corrected for both the host continuum (as mentioned 
in Sect. Hr) and Fe ii. 


Figure 3. Results from the [Fe ii] emission-subtraction. In the first piot, 
we .show the spectrum of J034740.i 8+010514.0 together with a fit of the 
[Fe ii] emission (yellow). In the second plot, shows the spectrum after 
suhtraction of the [Fe ii] emission. 



Figure 4. The plot shows the fit of the H/3 and [OIII] emission line 
complex with multiple Gauss functions for J0347. 


4.3 Emission-line fitting 

To ac curately measure th e fluxes of emission lines, we used 
MPFiTEXPR llMarkwardll 1 20091) . We fitted several Gaussians or a 
single Gaussian to the line profile, to fit narrow or broad compo¬ 
nents depending on the details of their shape or the presence of 
neighboring features. For example, to fit the Ha+iN ii] complex, 
we used four Gaussians components, where the flux ratio of the 
[N ii]/l 6583 A, 6548A doublet is fixed to the theoretical value 
of 2.96. In addition, their widths are supposed to be identical. 
Likewise, for the complexes of Yifi+[0 iii] A 5007 A, we fitted 
three Gaussians components as shown in Fig. [4] and two Gaus¬ 
sian components if the Hj3 line is narrow. Additionally, for some 
broad lines we fitted a Lorentzian profile if the width of the lines 
cannot be represented by Gaussian line profiles. The results of 
the line fitting are listed in Tab.|^ 


5 RESULTS 

Redshifts and luminosity distances Dy jHoggl [19991) are 
listed in Tab.[T] The luminosity distances Dy were calculated us¬ 
ing the redshift of the source combined with the cosmology con¬ 
stants assuming a Hubble constant Ng = 70 km s“' Mpc^', and 
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X offset [arcsec] x offset [arcsec] x offset [arcsec] 

Figure 2. The plot shows from left to right an SDSS image of the galaxy J1203, a star/point source from the same SDSS frame (coordinates are listed in 
Tab.|4) and the residuum that is left after subtracting the scaled star from the galaxy. The same plots for the other galaxies can be found in the appendix. 


Table 2. Comparison of seeing conditions for the sources for which we analyzed the spectra. 


Sources 

LBT seeing 

Comparing Spectro. 
Inst. seeing 

SDSS photometric 
Median Seeing (r-band) 

J1203 

0.82 

SDSS 

1.48 

1.4 

J0938 

1.30 

SDSS 

1.69 

1.4 

J0347 

1.40 

OHP 

2.5 

1.4 

J1158 

1.40 

SDSS 

1.22 

1.4 

J0802 

1.20 

SDSS 

1.49 

1.4 

J0911 

1.20 

SDSS 

1.23 

1.4 

J0354 

1.20 

Hiltner 

2.5 

- 

J0153 

1.90 

Beijing 

2.5 

- 


Table 3. The observed continuum flux density variability between LBT and SDSS/OHP of J0938, J1203, J1158, J0911, J0802, and J0347 from the 
photometry and spectroscopy aspect. ^"^For J0347 first epoch spectroscopy obtained by OHP. 


Sources 

Fill. 

Wa. 

A 

photometry 
SDSS 
[10“'* erg 
s“‘ cm“2A“'] 

Emission line 

contribution in 
spectral filters 
in % of continuum 

spectroscopy 
SDSS'“> 
[10“'* erg 
s“' cm“2A“'] 

spectroscopy 

LBT 

[10“'* erg 
s“' cm“2A“'] 


(u) 

3543 

1.26 ±0.03 


— 

2.04 ± 0.99 


(g) 

4770 

1.06 ±0.02 

9 

1.20 ±0.19 

1.54 ±0.38 

J1203 

(r) 

6231 

1.53 ±0.01 

25 

1.02 ±0.31 

0.87 ±0.13 


(i) 

7550 

1.09 ±0.01 

33 

0.81 ±0.21 

0.61 ±0.14 


(z) 

9134 

0.53 ± 0.03 


0.65 ± 0.32 

0.49 ±0.13 


(u) 

3543 

1.60 ±0.03 


— 

1.87 ±0.31 


(g) 

4600 

2.39 ± 0.02 

2 

3.41 ±0.19 

2.66 ± 0.27 

J0938 

(r) 

6231 

2.01 ± 0.02 


3.30 ±0.15 

2.48 ±0.31 


(i) 

7625 

1.67 ±0.02 

7 

3.05 ±0.17 

2.36 ±0.41 


(z) 

9134 

1.50 ±0.02 


2.79 ± 0.27 

2.12 ±0.45 


(u) 

3543 

68.81 ±0.01 


— 

52.79 ± 0.09 


(g) 

4770 

53.90 ± 0.02 

6 

44.60 ± 0.21 

60.38 ± 0.02 

J0347<“> 

(r) 

6231 

44.92 ± 0.01 

36 

31.12±0.15 

39.23 ± 0.01 


(i) 

7625 

40.91 ± 0.02 


27.71 ±0.11 

33.81 ±0.02 


(z) 

9134 

26.93 ± 0.01 


— 

26.39 ± 0.04 


(u) 

3543 

5.98 ± 0.04 


— 

5.12 ±0.63 


(g) 

4770 

4.29 ± 0.02 

5 

4.62 ± 0.23 

3.91 ±0.33 

J1158 

(r) 

6400 

2.87 ± 0.02 


2.81 ±0.14 

3.31 ±0.16 


(i) 

7700 

2.21 ± 0.02 

6 

2.15 ±0.16 

2.65 ± 0.03 


(z) 

9134 

1.77 ±0.03 


1.48 ±0.48 

1.87 ±0.07 


(u) 

3543 

0.41 ± 0.06 


— 

0.88 ± 0.49 


(g) 

4770 

0.97 ± 0.01 


1.02 ±0.12 

1.37±0.11 

J0802 

(r) 

6231 

1.17 ±0.01 

2 

1.14±0.11 

1.49 ±0.06 


(i) 

7625 

1.24 ±0.02 

3 

1.16 ±0.20 

1.71 ±0.10 


(z) 

9134 

1.04 ±0.02 


0.97 ± 0.25 

1.23 ±0.29 


(u) 

3543 

0.89 ± 0.04 


— 

0.61 ±0.27 


(g) 

4770 

0.86 ± 0.02 

5 

0.67 ± 0.04 

0.92 ± 0.22 

J0911 

(r) 

6200 

0.85 ± 0.01 


0.90 ± 0.04 

1.10±0.15 


(i) 

7700 

0.68 ± 0.02 

7 

0.72 ± 0.03 

0.96 ± 0.07 


(z) 

9134 

0.54 ± 0.03 


0.56 ±0.10 

0.76 ±0.19 
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Table 4. Compaiison of the apparent sizes of the sources and reference stars in the field. 


source 

Star 

RA Dec. 

(hh mm ss) 

FWHMl 

Star 

FWHM2 

Angle 

FWHMl 

Galaxy 

FWHM2 

Angle 

tr 

J1203 

12 02 59.4 

16 24 20.5 

0.95 ± 0.03 

1.03 ±0.03 

1.24 ±0.56 

1.26 ±0.15 

1.09 ±0.09 

10.97 ± 0.29 

J0938 

09 38 22.3 

13 49 54.4 

1.15 ±0.01 

1.04 ±0.01 

0.97 ±0.13 

1.31 ±0.07 

1.11 ±0.08 

1.88 ±0.48 

J0347 

03 47 35.8 

01 04 08.6 

0.99 ± 0.01 

1.31 ±0.01 

1.81 ±0.15 

1.09 ±0.01 

1.35 ±0.01 

0.87 ±0.16 

J1158 

11 58 14.2 

13 25 57.0 

1.43 ±0.02 

3.76 ± 1.48 

0.92 ± 0.30 

1.44 ±0.04 

1.49 ±0.04 

0.64 ± 0.48 

J0802 

08 01 38.2 

55 11 44.2 

0.57 ±0.01 

0.60 ± 0.01 

0.78 ±0.17 

0.67 ± 0.04 

0.52 ± 0.03 

49.44 ± 0.46 

J0911 

09 1! 47.7 

40 34 38.9 

0.75 ± 0.04 

0.81 ±0.04 

3.22 ± 0.82 

0.87 ± 0.06 

0.91 ±0.06 

0.79 ± 0.90 


Notes: For each galaxy, we give the coordinates of the star that we subtracted from the galaxy image as well as results of a Gaussian fit to galaxy and 
star. We used an elliptically shaped Gaussian function. 


Table 5. Continuum fits. 


name 

instrument 

Ay (stars) 

stars 

powerlaw 

Fe II 

J1203 

LBT 

0 

60% 

40% 

0% 


SDSS 

0 

60% 

40% 

0% 

J0938 

LBT 

1.5 

90% 

10% 

0% 


SDSS 

1.5 

80% 

20% 

0% 

J0347 

LBT 

1.2 

57% 

20% 

23% 


OHP 

0.8 

72% 

13% 

15% 

J1158 

LBT 

0 

0% 

85% 

15% 


SDSS 

0 

0% 

85% 

15% 

J0802 

LBT 

0 

50% 

30% 

20% 


SDSS 

0 

52% 

33% 

15% 

J0911 

LBT 

0 

50% 

40% 

10% 


SDSS 

0 

45% 

45% 

10% 

J0354 

LBT 

1.8 

60% 

30% 

10% 

JO 153 

LBT 

0 

62% 

23% 

15% 


Notes: For each observation we give the flux contributions of the stellar component, the AGN/powerlaw component, and the Fe ii template, as well as 
the extinction Ay of the stellar component in the wavelength interval 5100 A to 5600 A (restframe wavelength). 


a standard cosmolog y with parameters = 0.3, and Ha = 0.7 
iSoergel et alJl20o3 !^ which we use throughout the paper). 

The observed continuum flux density variability of the sources 
derived from LBT and SDSS data is summarized in Tab. (3] 

The main body of comparative variability studies of extra- 
galactic nuclei is carried out without a detailed dilferentiation 
of contributions of different spectral components to the overall 
emission. For the brighter QSOs this can be justified since the 
host contribution is small. Such a spectral analysis requires a 
very high signal to noise in order to determine the exact nature 
of the stellar host contribution which in return is dependant on 
the modeling details of the dominant contributing stellar popula¬ 
tions. Given the inter-calibration quality for the newly presented 
8 sources such a analysis shows that around 5500 A (restframe 
wavelength) the stellar contribution to the integrated flux for the 
S2/NLS1 source is >50% and for the QSOs it is ^50%. How¬ 
ever, the quality of the data does not allow to study the variabil¬ 
ity of the power spectrum component alone, since the details of 
the stellar population analysis required a much higher data qual¬ 
ity. Hence, the difference in the continuum variability properties 
between S2/NLS1 and QSO/BLS1 presented in this paper is in¬ 
fluenced by the presence of a stronger stellar contribution for the 
S2/NLS1. Details are given in column 5 of Tab. [3] If no value is 
given in column 5 then the contamination is below 5%. 

The variability of a larger sample (see section IbTt is dis¬ 
cussed with continuum and line variability extracted from spec¬ 
tra. In order to avoid contamination from the host and the Fe¬ 
lines, the line variability is measured from spectra corrected for 
these contributions (see subsections above). However, in order 
to be able to enlarge the sample with data from the literature, the 


continuum variability measurements do contain the stellar flux 
contribution. The effects of this are discussed in section [63l 

The line fitting results are listed in Tab. |9] We plot all our 
spectra in the rest wavelength as shown in Fig. i5] |A^lA9l etc.). 
We did that by applying /!„„ = /toj,.,/! -l- z. Another correction 
was applied to the flux density (cosmological dilution), changing 
this value to the rest system via 

Fa rest = Fa obs. x (1 -l- z)^ . (1) 

All continuum flux densities and luminosities listed in the tables 
are derived as observed quantities (i.e. without corrections). The 
listed line fluxes and luminosities have been corrected for Fell 
and stellar continuum as described above. 

5.1 SDSS J093801.63+135317.0 

SDSS J093801.63-tl35317.0 (hereafter J0938) is discussed 
in the discovery of a population of normal field galaxies that 
have luminous and strong FHI L (forbidden high by ionize d 
lines) and HeII/14686 emission dStenanian & Afanas’ evil 201 ih . 
The authors used the ’Bolshoi Teleskop Alt-azimutalnyi BTA- 
6' 6-m diameter teles cope in Russia and obtained s pectra with 
0.86A/px resolution. IStenanian & Afanas ’evl ll201lh report, in 
their Table 2, the value of FWHM and flux for emission lines in 
the optical spectrum of J0938 taken with two instruments with 
the same spectral resolution ofR~ 2000 at two different epochs. 
The first is from the SDSS telescope on 18.12.2006 and the sec¬ 
ond is from the BTA-6 on 16.04.2010. 

lYang et ^ ll2013h present results of seven rare extreme 
coronal line emitting galaxies, and 10938 is also among these 
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sources. These objects are reported by IWang et ^ ( l2012h . and 
four of these galaxies (J0938 and three others) have a large 
variability in coronal line flux, making them good candidates 
for tidal disruption events (TDEs). They detected a broad He 
11/1468 6 emission line in the spectrum of J0938. lYang et aP 
(l201.3h find broad coronal and high-ionization lines that are su¬ 
perimposed on narrow low-ionization lines. They interpret this 
finding as indication that J0938 is a composite of a S2 nu¬ 
cleus and a star formation region. In the SIMBAD catalogue 
the object is listed as an H ii galaxy. However, according to the 
[O i]/i6300/Ha versus [O iiil/i500 7/Hg diagnostic diagram in 
IWang et ^ (1201 2h and I Yang et al.l (l201.3h it’s appearance to be 
that the source J0938 is located in the region between H ii and 
S2 galaxies (here we adopt an HII/S2 composite nucleus clas¬ 
sification) and belongs to a sample of sources with spectra that 
are dominated by the interaction of a super massive black hole 
with the nuclear environment l Stenanian & Afanas ^ l201ll: 


IWang et alJl2012l:lYang et alj201^ 

The variability data for this source are consistent with the 
uncertain classification of the source as HII/S2 and support the 
presence of a variable nucleus. In Tab.l^we list the flux density 
measurements of J0938 for three epochs at different wavelengths 
as obtained through our LBT observations and via the spec¬ 
troscopy and photometry results listed by SDSS. We note that 
flux values of the continuum of LBT spectroscopy are in good 
agreement with photometric data of SDSS at the same wave¬ 
lengths, while both are different to what can be derived from 
SDSS spectroscopy. We notice that their are of order 30% vari¬ 
ations in line and continuum flux between the different epochs 
we use for our investigation. 

As we show in Fig. lAll and as is evident from the data in 
Tab.|4] the host of the source is extended. Comparison with stars 
in the field show that about 50% to 60% of the continuum flux 
in the LBT 0.8” slit is due to the unresolved nuclear component. 
The rest can be attributed to the a contribution of the extended 
host. Hence the estimated continuum variability is probably only 
an upper limit. However, for completeness we leave the variabil¬ 
ity estimates in Fig. as upper limits. The median and median 
deviation derived from this plot are not effected by this. In Fig.[5] 
we over-plot the two optical spectra (LBT and SDSS) of J0938. 
The LBT spectrum was observed in 2012, and the SDSS spec¬ 
trum was taken in 2006. Therefore, we can discuss the variabil¬ 
ity on a time scale of 6 years. We found that, for this object, 
there are two types of variability: first, there is variability in the 
continuum level between the two spectra, where the continuum 
level of the SDSS spectrum is higher by a factor of ~L3 than 
that of the LBT spectrum. This can be seen in Fig.[5](a,b,c,and 
d) where we plot different sections of the J0938 spectrum. The 
second type of the variability is in the emission lines. To show 
this we subtract the two spectra from each other and plot the dif¬ 
ference (Fig.|5]e, f, and g). Additionally, we show the variability 
of these lines by displaying the ratio of the two spectra (Fig. [5] 
h,i, and j). We found that the emissions lines that varied most 
in the spectrum of J0938 are [O ii], H/3, [O iii]/i5007, and Ho-. 
The Ha and [O iii] lines show a variation of the order 1.25. The 
measured line fluxes are listed in Tab.|3 


5.2 SDSS J120300.19+162443.8 

SDSS J 1203Q0.19-hl62443.8 (hencefo rth J1203) has been 
mentioned in IShirazi & BrinchmannI ( l2012h among 2865 galax¬ 
ies to have a strong nebular He 11/14686 emission. A strong 
He 11/14686 line indicates that the nuclear radiation field of these 
objects is dominated by highly ionizing radiation. The ioniza¬ 
tion potential of He^ is 54.4 eV corresponding to a UV photon 
wavelength of /I a; 228A. 


In Tab. [3] we list the flux densities of J1203 at different 
wavelengths obtained from our LBT data and spectroscopy and 
photometry data as listed by SDSS. We note that the flux values 
of the continuum at different wavelength indicate variability. 
Moreover, the power law index of the continuum feature of 
both observations SDSS (photometric & spectroscopy) and 
LBT (spectroscopy) has also varied as we show in the optical 
spectrum of J1203 in Fig. lA3l 

The optical emission line spectrum of J1203 is dominated 
by a NLR (see Fig. IA3t . In Tab. we see that all the emis¬ 
sion lines are narrow with FWHM values of < 8.3A (where 
8.3A s! 500 km s“'). J1203 has very low count rates at the 
continuum level. The Ha c omplex gives indica tions of a possi¬ 
ble faint broad component. iBarth et ^ (l2014h find that all S2 
candidates in their sample for which high spectral resolution 
deep exposures had been taken showed some broad line emis¬ 
sion that was not visible in previous spectra. A comparison be¬ 
tween the [OIII] /I5007 and H/f line shows a very similar line 
profile, within the uncertainties, hence we classify J1203 as a S2 
galaxy. 

The over-plot of the two optical spectra (LBT and SDSS) 
for J1203 is presented in Fie. lA3l The SDSS spectrum was ob¬ 
served in 2007, while the LBT spectrum was observed in 2012. 
Hence, for this galaxy we can discuss the variability that hap¬ 
pened over 5 years. We find that the continuum variability here is 
rather limited to the blue channel region from 3300A to 3700A 
where effects o f red wing of t he ’’Big Blue Bump BBB” ma y 
well be present (lGaskellll2008l : IStarling & Puchnarewiczll200lh . 
as shown in Fig. IA3l fd). The observed continuum of the LBT 
spectrum in this region is higher than the SDSS spectrum. But 
for the rest of the spectrum both spectra are at the same level till 
6000 A. The line emission varies by about 40% (see Tab.|9). 


5.3 SDSS J115816.72-h132624.1 


The optical spectrum of SDSS J115816.72-tl32624.1 
(hereafter J1158; observing dates see Tab. [T) is strongly dom¬ 
inated b y Fe ii lines (as can be seen in Fig. |A4lan d as reported 


also in 


Kovacevic. Ponovic & Dimitriievid boid : IPong et ^ 


l201lh . Tab. [3] shows that the flux values of the photometric 
and spectroscopy observations obtained within the SDSS survey 
are in good agreement. However, both values show a difference 
with respect to the flux v alue extracted from the LBT spect rum, 
indicating variability. In IWang & Rowan-RobinsonI ( l2009lf the 
authors refer to JII58 as a faint source and mention it in the 
Imperial IRAS-FSC Redshift Catalogue (IIFSCz) among 60303 
galaxies. Verifying the consistency of source positions they im- 
prove d the optical , near- infrared, radio identification. Accord¬ 
ing to lZhang et al.l ll2010h J1158 is classified as a low-ionization 
broad absorption-line (loBAL) quasar, with [Mg ii] absorption¬ 
lines with a width of Avc > 1600 km s“*. Additionally, the au¬ 
thors calculate the continuum luminosity at 5100A (/iLsioo) as 
7.7 X 10‘'''erg s“‘. From the LBT spectrum we obtain a luminos¬ 
ity at the same wavelength with a flux of 2.81 x 10"^^erg s“'. The 
SDSS spectrum results in a value of 2.25 x lO'^^erg s“*. Con¬ 
sistent with the fluxes listed in Tab. [3] we can conclude that the 
continuum luminosity of this source has varied by a factor of 1.2 
over the past 8 years. 

In Fig. lA4l we present the over-plot of the LBT and SDSS 
spectra for J1158. The LBT data have been obtained in January 
2012, and the SDSS survey monitored this source in May 2006. 
Hence, for this object we are able to study and discuss the spec¬ 
tral variability on a time scale of 6 years. After over-plotting 
both spectra and measuring the flux densities of the emission 
lines (Tab. |9) for this galaxy, we found that the variability in 
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Figure 5. The plot shows the optical spectrum of J0938 as observed by SDSS and LBT-MODS. The second row shows zooms into the spectrum in 
different regions. The third row shows the differences between SDSS and LBT spectrum in these regions, while the fourth row shows the ratios of these 
spectra. 
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the continuum level and the continuum spectral index (hereafter 
PLI) may well be coupled to each other as the red tail of the BBB 
may be the dominant contributor to the variability behavior. In 
addition we find strong variability in emission line intensity. 

As shown in Fig. lA4l (a.b.c. and d) the continuum level of 
the LBT spectrum is higher than that in the SDSS spectrum by 
a factor of about 1.25. This could be linked to power-law vari¬ 
ability over the entire optical spectrum as expected from local- 
ized temperature fluctuations of a simple inhomogeneous disk 
iRuan et alj2014ll . We find the PLI for both spectra to vary from 
-0.21 erg s“' cm“^ / A for the SDSS spectrum to the PLI of LBT 
spectrum of -O.OI erg s“* cm“^ / A. Among them the H6/I486I 
and [O iii]/l5007 line show the strongest variation by a factor of 
about 3. In general, the emission lines in this source are more 
variable than in other objects of our sample. 


metric and LBT spectroscopy are similar they both are different 
to what can be derived from SDSS spectroscopy. This indicates 
that J0802 shows some continuum variability between epochs. 

The LBT spectrum was observed in January 2012, while 
the SDSS spectrum was taken in January 2005, resulting in a 
time baseline of 7 years for this source. Plotting both spectra 
(SDSS & LBT) for J0802 in Fig.|^(a, b, c, and d) we find that 
there are spectral similarities to the source J0938. Both spectra 
have the same continuum slope but the continuum level of the 
LBT spectrum is higher by a factor of 1.3 compared to the SDSS 
spectrum. The intensity of the emission lines in LBT spectrum 
are stronger than those obtained from the SDSS spectrum (see 
Tab. 13 • The Hy 44340 line shows the strongest variation with a 
factor of about 2. 


5.4 SDSS 1091146.06-1-403501.0 

SDSS J091146.06-H403501.0 (hereafter J0911) was first 
reported in iMcMahon et al.l l l2002h as a radio source. In the 
third data release of the S PSS survey l lSchneider et al.ll200.5h 
it is classified as a quasar. IZhang et al.l | |201(]|) classiW J0911 
as a low-ionization broad absorption-line (loBAL) quasar, with 
[Mg ii] absorption-line in its spectrum with a line width of 
Auc > 1600 km s^L J0911 exhibits strong Fe ii lines as 
shown in Fig. IA5I Based on the SDSS DR5 data IZhang et alJ 
calculate the continuum luminosity of J0911 at 5100A 
(/ILsioo) to be 2.9 x 10'*'*erg s“’, while at the same wavelength 
we derive from our LBT spectroscopy data a luminosity value 
of 3.66 X 10''‘*erg s“*, and for the SDSS spectrum we find 
2.84 X 10‘'‘*erg s“*. These three values, and the measurements 
in Tab. [3 show that there is significant variability in the contin¬ 
uum emission of this source over these two epochs (see Fig. IMJ- 
This is in contrast to the strong variability we find for J1158 and 
supports the consistency of the calibration between the different 
data sets. As part of the SDSS survey this source has been ob¬ 
served in March 2003, while the LBT spectrum of the source 
was taken in February 2012. Hence we look at a time span of 9 
years for the variability search. The optical spectra of J0911 can 
be seen in Fig. lA5l (a. b, c, and d). As illustrated in this figure, the 
variability of the source is small compared to other sources in the 
sample, except for the blue part of the spectrum in Fig. lA5l (b). 
where we see a difference in the continuum slope. This variation 
in the blue part of the continuum may caused by an accelerat¬ 
ing outflow emanating from the blac k hole in the center of the 
galaxy, see lShapovalova et al.l l l2010t) . In Tab. we list the flux 
measurements of the emission lines and their FWHM. 


5.5 SDSS 1080248.18-1-551328.9 

SDSS J080248.19-H551328.9 (in the following J0802) was 
first mentioned in the Fifth Data Rele ase DR5 of the SDSS cata¬ 
log for quasars iSchneider et alJ2007h . T he initial redshift deter - 
mination z = 0.66287 ± 0.00107 in DR5 jSchneider et alj2007h . 
later has been correcte d to a value of z = 0.664 065 ± 0.000355 
He^tt & Wildll2oI(ih . I Schneider et alJ i2010h and llnada et ^ 
l2012t) class i fied J0 8 02 as a quasar, which is su pported by 
iGibson eTa I] l l2009ll . iLundgren et alj l l2009ll and I Allen et Si 
Imil) with reference to a broad absorption line of [Mg ii] 42800 
in its spectrum. From the optical spectra it is apparent that J0802 
is highly dominated by Fe ii lines (see Fig. lA6b . We notice that 
the continuum power law index derived from the LBT and SDSS 
spectroscopy data as well as the value derived from SDSS pho¬ 
tometry are in good agreement (Tab.[3j- According to these data 
the continuum spectrum peaks in the i-band. However, compar¬ 
ing the flux densities we find that while values of SDSS photo- 


5.6 2MASX J035409.48H-024930.7 


This source 2MASX J035409.48-f024930.7 (hereafter 
J0354) is one of the two sources that have not been cov¬ 
ered by the SDSS survey. J0354 was serendipitously dis- 
covered and discussed, to g ether with 18 other objects, in 
IChanan ^^PovTOS & Margotil ( Il98lh as an AGN X-ray source. 
In Bothun et al. lll982li the authors report that the J0354 spec¬ 
trum (Mv ~ 18) allows for either QSO or Seyfert 1 classi- 
fication. They outline that comp ared to other samples (e.g., 
iMargon. Chanan & DownesI Il982h the object is amongst the 
brightest X-ray sources. They also report that this source has 
a large ratio of X-ray to optical luminosity Lj,/L„p, = 2 
which is above the t y pical ratio of 0.5 found for quasars. 
iHaddad & VanderriestI il99lh state that J0354 is a gas rich 
galaxy that arose from a two-galaxy interaction, where the larger 
one is either a quasar or SI and the second clearly shows 
signs of tidal distortion. iHutchings. Campbell & Cramptonl 
198^ suggest that the com panion is bluer than the quasar. 


IVeron-Cettv & Veronl(l2006li classify the source as SI.5. 

We can compare our LBT spectroscopy data for J0354 with 
spectroscopy data obtained in Oct 2004 by the 2.4 m Hilt- 
ner Telescope at MDM Observatory at Ki t t Peak , Arizona, 
USA published by lOrupe. Pradhan & Frank! (l200.5h . Here the 
source was observed with 1.5” and 2” slits under moderate 
seeing conditions. We notice that the continuum level of the 
LBT spectrum is different than the spectrum of Hiltner Tele¬ 
scope as shown in Fig. lA7l where the continuum level of J0354 
taken with the Hiltner telescope is higher than the LBT spec¬ 
trum for the same source, and the emission lines of their ob- 
servation are more inten s e than our observation. Additionally, 
iGrupe. Pradhan & FrankI ll200.5h report line fluxes of [[O iii] 
44959 and 45007 as (130±6) and (413± 15)x lO^'^erg s^' cm^^, 
respectively. From the LBT data we obtain fluxes of the same 
lines as (90) x 10“'*erg s“' cm“^ and (505) x lO^'^erg s“' cm“^, 
respectively. The comparison shows that the [O iii] 45007 line 
varies by a factor of almost 50%. 


5.7 GALEXASC J015328.23+260938.5 

GALEXASC J015328.23-H260938.5 (henceforward J0153) 
is the second source that is not covered by the SDSS sur¬ 
vey. This source is present in the I nfrared Astronomy S atel- 
lite (IRAS) Point Source Catalogue llBeichman et alJfl988l . us¬ 
ing the name IRAS 015064-2554) and in the 2nd XMM- 
Newton Serendipitous Source Cat alog jMemola et alj l2007l : 
iNoguchi. Terashima & Awal3 l2009l using the name 2XMM 
JO 15328.44-26093). From the IRAS survey, we know the flux 
density at two wavelengths, 60 and 100 microns, to be 0.5777 & 
1.144 Jy, respectively. In addition to our LBT obser vations (see 
Fig. lA8t a spectrum of J0153 has been obtained bv IZhou et alj 
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Table 6. Median and median deviations for the variability in the contin¬ 
uum & emission lines. 


Type 

Continuum 

NaiTow 




Emission lines 


p 

dfi 

p 

dfi 

BLSl&QSO 

39.0 

17.0 

30.0 

10.5 

S2&NLS1 

27.0 

9.0 

15.0 

6.0 


Notes: fi are the median values, dfi the median deviations as defined in 
the text. 


Table 8. The vaiiability percentage in continuum at rest wavelength 
5100Aand Wfi, for AGNs taken from reverberation mapping (light 
curve). 


Sources 

Cont. 5100A 

% 

% 

Type 

References 

Mrk 335 

34 

24 

NLSl 

1,2 

Mrk6 

36 

25 

NLSl 

1 

3C 120 

57 

40 

BLSl 

1 

Ark 120 

48 

47 

BLSl 

1 

Mrk 590 

43 

50 

BLSl 

1 


References. (1 ) \Peterson et al\ \l99A) . (2} \Grier etal.\ \20li} . 


120021) in February 2002 with the Zeiss universal spectrograph 
located at the Beijing Observatory 2.16 m telescope, with a slit 
width of 2'.'5 and seeing disk of 2'.'5. Moreover, the spectral res¬ 
olution determined on the night sky is 5.2A FWHM. They also 
found that J0153 is a radio-loud quasar with very strong Fe ii 
emission lines, for the far- and near-infrared luminosity, they 
found LpiR » 10*^'^ Lq and Lnir a; 10‘^'^ Lf?). 

JO 153 is one of 30 galaxies that ICombes et alJ 1201 ih study to 
follow the galaxy evolution and especially the star formation ef¬ 
ficiency (SFE) using CO lines. They find that the SFE of this 
source is > 802 (Lq/Mq) which is 4.7 times higher than the 
local ultra-luminous infrared galaxies (ULIRGs) (170 Lq/Mq), 
indicating highly efficient star formation activity. 


5.8 SDSS J034740.18-H010514.0 


SDSS J034740.18-H010514.0 (henceforth J0347) was men- 
tioned first tim e in a study for high-luminosity sources 
ILow et alJI 19881 using the name IRAS 034504-0055) and clas¬ 
sified as quasar. In|HewitL&&urbidgd 11991 ) J03 47 was classi¬ 
fied as S1 ga laxy bv1vdroi>CettT^ V6rorJT 2006l) as SI.5. With 


10‘°^ lLowetal.lll988h and| 


Low et al 


1 19^ find .10347 


Lfir : 

to have the lowest luminosity among other quasars in their sam¬ 
ple, additionally they find that this source has an extremely red 
optical continuum. Since no SDSS spectroscopic observations 
are available for this object, we used spectroscopic data taken on 
21 September 2003 with the CARELEC spectrograph attached 
to the L93m telescope of the Observatoire de Haute-Provence 
(OFIP) with a 2" slit under 2.5" seeing. 

iGiannuzzo & Stir^ i ll9961) search for variability in a sample of 
12 narrow line Seyfert 1 galaxies, including J0347. They find 10 
of these sources to be variable in the flux of the optical permit¬ 
ted lines over a period of one year and also found variation in the 
continuum level. The variability ratio they find in the permitted 
lines and Flo' is 4.58 and 0.599 respectively. 

We list measurements of the continuum flux from SDSS pho¬ 
tometry as well as from OFIP and LBT spectroscopy in Tab. [3 
Line flux measurements are presented in Tab.|3(see also plots in 
Fig.lA9t. 


6 DISCUSSION 


Our study allows us to compare the continuum and line 
variability of two different source samples with different spec¬ 
troscopic identifications that by themselves imply a different de¬ 
gree of activity. We find that the Seyfert 2/NLS 1 sample shows 
in general a smaller degree of variability compared to the BLS1 
and QSO sample. In the following we probe the significance of 
this difference and present a possible physical interpretation of 
the observed phenomenon. We extend the data on the contin¬ 
uum variability with multi epoch i nformation on ESQ 0 1 2-G2L 
NGC 1365, and MKN 1044 from IGiannuzzo & Stirnd lll99fil) . 
NGC 7603 from iKollatschnv, BischoT^^ietricinTOOOl) and 


PKS 2349-14 from IKollatschnv. Zetzl & Dietrichl ( l2006al) . We 
extended the data on multi epoch line varia bility with ESQ 012- 
G2L NGC 1365 and MKN 1044 from IGiannuzzo & Stirnel 
l ll996l) . 


6.1 Probing the degree of variability 


Based on our detailed study of 18 sources (see Tab.[T] Tab.[8] 
and captions in Figs. and [7} we find that there are signif¬ 
icant differences in variability between S2(&NLS1 (8) sources 
and BLSli&QSO (10) sources. This applies both for the con¬ 
tinuum and the emission lines. In the following we discuss the 
relative variability A„,ov as the maximum difference 6„,a:, with 
respect to the maximum value S„,ax, i-e. A„,aj- = 6„axlSmax in 
(i.e. dependent on the inter-calibration uncertainties which are 
the same for both source classes). Correction for relative mea¬ 
surement uncertainties <t of the order of 10% results in 


A “ /a^ — rr^ 

^ma.x,corr. yj ‘-^max ^ 


( 2 ) 


We find t hat the BLSl&Q SO have stronger variablity 
than S2&NLSL lAi et akl ( l2013ah came to the same conclusion: 
NLS1 have a systematically lower degree of variability if com¬ 
pared to BLSl sources. This is in agreement with the existing 
anti-correlation between AGNs varia b ility and Eddingt on ratio 
iMeusinger. F inze^^ Hooni (I 2 OIII) IZuo et alJ 1 I 2 OI 2 I) . Other 
authors (e.g., ^^bb'&^aikMDOOOl) concluded from variability 
studies that the blue part of AGN spectra is relatively stronger 
variable than the red part. For our study we show this result in 
Fig. [6] We notice that for J0354 and 10911 the variability in the 
blue part of the continuum is, larger than the variability in the 
red part. 

To judge on the difference between the two source samples 
in their continuum and line variability we calculated the median 
and median deviation as well as the mean and standard devia¬ 
tion. We also chose the median since we do not know ab-initio 
if we have normal distributions. Here, the median deviation is 
defined as the median of the absolute differences between the 
median and the sample values. The median deviation is usually 
smaller than the standard deviation as it is less sensitive with 
respect to outliers. From Tab. and Figs. and [7] we find that 
the median degree of continuum variability is 27%±9% for the 
S2/NLS1 and 39%±17% for the BLSl/QSO sources. For the 
median line variability we find 15%±6% for the S2/NLS1 and 
30%± 11 % for the BLS 1/QSO sources. Furthermore, we find that 
the median values are separated by about one median deviation. 


6.1.1 Robustness of the separation 

Results for narrow and broad line components: In Fig. [7] 
we show the results for the combined narrow and broad line 
recombination line fluxes in order to compare the data with 
literature values. However, the difference in the variability 
behavior between the narrow line and the broad line sources 
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Table 7. The T-test statistic for the variability in the continuum & emission lines. 


Type 



Continuum 




Emission lines 



n 

T 

(T df 

T-value 

n 

T 

(T 

df 

T-value 

BLSl&QSO 

42 

42.25 

20.76 

11.72 

3.38 

18 

29.48 

10.71 

33 

4.25 

S2&NLS1 

35 

29.51 

17 

15.51 

8.67 


Notes: Here n denotes the number of sample elements, x ^t-e the sample mean values, cr the sample standard deviations, df the degrees of freedom 
followed by the corresponding T-values. 
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Figure 6. The variability percentage in the continuum. The data for the NLSl sources ESO 012-G21, NGC 1365, and MKN 1044 are taken from 
Giannuzzo & Stirpe (1996). The data for the BLSl sources NGC 7603 and PKS 2349-14 are based on spectra from Kollatschny et al. (2000) and 
Kollatschny et al. (2006a), respectively. 


(with 15%±6% and 30%±11%, respectively) also holds if 
we investigate the narrow line and broad line components we 
obtained from fitting the hydrogen recombination lines (see 
Tab.lAB. Here the median value for the variability of the narrow 
line objects J0347, J1203 and J0938 of 20±4 lies well below 
the median result for the broad line objects of 41±8 derived 
from the 4 objects J1158, J0911, J0002 and J0153. Although the 
separation into the two different source classes has been done 
using an overall line width criterion, we find that the difference 
is also present if we look at the NLR and BLR line components 
separately (see discussion in section l63l l. 

Statistical considerations: For the BLSl/QSO sample we 
see the tendency in some sources for the variability to be up to 
20% stronger in the blue (3000A) compared to the red (close 
to 6000A). However, this tendency is not equally well fullfiled 
for all sources and rather weak with respect to the overall 
variability range from about 20% to 70%. However, here we are 
more interested in the overall continuum and line variability. 
Both are estimated for different sources and in different lines 
or spectral regions of the continuum, which are all subject to 
different wavelength and line strength dependant calibration 
uncertainties. Hence, we assume that for the two source samples 
the variability estimates are sufficiently independent, and can 
be represented by a mean distribut i on. In this case we can 
apply a T-test dSenn & RichardsonI 1 19941) to determine the 
probability that the data sets from both the continuum and line 


variation originate from different distributions. The relevant test 
quantities are given in Tab.|7] 

We applied a T-test for two sample experimental statistics for 
independent groups. The T-value was calculated via 

T = . (3) 

crdiff 

The standard error of the difference we calculated via 

Since we do not know the population standard deviation we 
use the sample standard deviation to estimate the standard error: 



n being the number of sample elements. The number of de¬ 
grees of freedom df as used for two independent groups we ob¬ 
tained as 

df = (m - 1) -H («2 - 1). (6) 

The probabilities that describe the matchi ng of the two 
samples are tabulated dSenn & Richards^ 1 19941) as a function 
of the T-values and the degrees of freedom df. We find that the 


© 0000 RAS, MNRAS 000,[T]Q7] 







































































Variability in Active Galaxies 15 


100 


to 

OJ 

c 

c 

o 


E 

OJ 

OJ 


o 

OJ 

fD 

c 

OJ 

u 

OJ 

Q. 




80 


60 - 


40 - 


~ 20 


Narrow line Objects. 


itiedian'deViatioh' 


• J1203 
A J0g83 
t J0347 

* ESO 012-G21 
■ NGC 1365 

O MKN 1044 


- 1 --^ ' 

:• 1 ^ 

• 


1 

1 

1 

1 

■ 



u 

i 

t 

i 












Broad. Line. Objects.. 


median deviation 


• J0354 

• J0802 

■ jogii 

• J1158 

• J0153 








cC'- 






sO'' 


r# 






N<^ 










9^ 








Figure 7. The variability percentage in the emission lines. The line measurements of ESO 012-G21, NGC 1365, and MKN 1044 are taken from 
Giannuzzo & Stirpe (1996). 


continuum and line variability of the BLSl and QSO sample 

originated with a better than 99.9% probability from a different ^ 

distribution compared to the Seyfert 2 and NLSl sample. Lcom. M x—M . (7) 


Comparison to other samples: The separation between the 
narrow line and the broad line sources is best fulfilled for the 
line variation. As shown in Tab. [8] our results are very much 
consistent with those for a set of NLSl and BLSl sources 
obtained from reve rber ation measuremen ts carried out by 
|Peterson~ Ii] h998t) an d I Grier et alj ( l2012h . 

The findings for our small sample are also in full agreement 
with the results of a comparative study of the optical/ul t raviole t 
variability of NLSl and BLS 1-type sources bv lAi etaL|j2013bh 
(based on 55 NLSl- and 108 BLSl-type nuclei). This underlines 
the representative character of our sample. The authors show 
that the majority of NLS 1-type objects exhibit significant 
variability on timescales from about 10 days up to a few 
years, but on average their variability amplitudes are much 
smaller compared to B LSl-type sources (see also, lYin et al] 
I 2 OO 9 I : iBarthet alJl2014h . Hence, in summary, we assume that 
the difference in continuum and line variability between the 
two samples is sufficiently significant to search for a physical 
explanation for the phenomenon. 


6.2 Accretion dominated variability 

Here we investigate how the difference in narrow line vari¬ 
ability between the objects with spectra dominated by narrow 
lines and broad lines is linked to the continuum variations of the 
active nuclei in the different samples. 

In the following we denote with A differentiation with re¬ 
spect to the different samples, i.e. A = , . If the variability 

is due to mass accretion onto a super-massive black hole 
with mass M we can write the continuum luminosity as 


The expected difference in continuum variability between 
the samples can then be expressed as 

d d 

AL„,„, oc AM X —M + M X A — M . (8) 

dt dt 

For recombination lines the variation of the nuclear contin¬ 
uum results in a variation of the number of atomic species that 
reach the next higher state of ionization out of which they can re¬ 
combine. For collisionally excited states this next higher state of 
ionization is the state in which they can be observed at low den¬ 
sities emitting forbidden narrow line radiation. In the following 
we use the case of the recombination lines to discuss the effect of 
the reverberation. Since we are mainly interested in the variation 
of the narrow line luminosity L/,„e as a reverberation response to 
the variations in continuum luminosity we can write: 


Cline ^ 


hVline Oeff d 
(hv) ttg dt 


M = e—M-, 
dt 


(9) 


Here vune is the line frequency, Ueff and ag are the effec¬ 
tive and MenzeFs case B specific recombination coefficients and 
(hv) denotes the mean energy per photon. Through the latter the 
combined coefficient ^ depends on the ration between the total 
illuminated cross-section Ac and ionized volume Vc of the re¬ 
verberating clouds in a single galactic nucleus. This ratio results 
in a size I of the reverberating volume which has a square root 
dependency on the continuum luminosity: 


^ OC ^ ocloc ^/Lcont. ■ (10) 

Therefore, combining the above equations for the line lu¬ 
minosity we can write: 
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L,,„, cx . (11) 

at 

The expected difference in line variability between the sam¬ 
ples can then be expressed as 



We now assume that - as a second order derivative - the 
variation A^M of the accretion stream onto the super-massive 
black holes from sample to sample is small and a finite difference 
AM in total mean black hole mass per sample needs to be con¬ 
sidered. In such a scenario, the strength of the accretion stream 
as such is the dominant quantity for the variations in continuum 
and line luminosity between between different sample members 
and the two sources sample in general. We can then neglect the 
second terms in equations and and write for the difference in 
continuum and line luminosity: 


d d 

AL,„,„ oc AMx—Moc —M (13) 

dt dt 

ALu„e « AM{^M)i oc (^M)i . (14) 

at at 

Hence the increase in activity in continuum and line emis¬ 
sion is due to the difference in AM between the samples, how¬ 
ever, there is an additional modulation for the variability of the 
line emission that is a function of ^M. Combining these results 
we find: 


^Eine ^ (^Eon!.)^ • ( 15 ) 

This simple relation was obtained by assuming that the dif¬ 
ference in accretion streams between the two samples is - to first 
order - negligible. The accretion stream onto the super-massive 
black hole is probably closely tied to the properties of the inter¬ 
stellar matter in its immediate vicinity. Hence, it may be different 
from source to source and may not strongly depend on the mass 
of the central super-massive black hole. However, if A^M needs 
to be considered it will have an influence on the general trend of 
the derived relation (i.e the overall slope and curvature of the 
trend). In fact, it may be responsible for part of the scatter of 
data points about this relation. However, we assume the second 
order derivative to be of lesser importance than the first order 
derivatives. This is supported by the fact that the general trend 
of the distribution of data points in Fig.[8]is well described. 

Looking at the median and mean values for the different 
samples in Tabs.|^and[7]we find that this condition described by 
equation [15] is to first order fullfiled to within about 15% : For 

3 

the means of the continuum variability we find (39/30)i = 1.5 
whereas the measured ratio of the means of the line variability 
is 27/15 = 1.8. For the median values of the continuum vari- 

3 

ability we find (42/35)i = 1.3 whereas the measured ratio of 
the median values of the line variability is 29/16 = 1.8. This 
proportionality relation between the continuum and narrow line 
variability is demonstrated in Fig.[8]and discussed in the follow¬ 
ing section. 


6.3 NLR response to the nuclear continuum variability 

In Fig.[8]the black crosses represent the median values and 
their uncertainties as shown in Figs.[^&[71 In this plot the curved 
lines represent the relation AL;,,,^ = e(AL„„,.)i, with e=1.0 for 
the thick straight line which follows our data surprisingly well. 
There are two main contributers and to that value such 


that e = For our sample, of the order of 50% of the con- 
tinuum flux density is likely to be due to the stellar contin¬ 
uum, which is not variable on the time scales discussed here. 
For more luminous objects this may be of the order of 20% 
dKollatschnv. Zetzl & Dietrichll2006a) . This indicates that the 
pure nuclear continuum variability which is not contaminated 
by stellar flux may in general be larger by up to 50%. This re¬ 
sults in a value of 77^,^^ ~1.9 (dashed line in Fig. [8). However, 
the covering factor rjfm of the clouds that are illuminated by the 
nucleus can be of the same order of magnitude. Hence, with a 
considerable scatter, e may in fact be close to unity, such that 
the bold black line in Fig.[8]that matches the observed data is an 
acceptable representation of the relation between the continuum 
and narrow line variability. 

We have carried out the interpretation of the line re¬ 
sponse using the reverberation formalism. This is usually ap¬ 
plied to broad lines only and used to derive the size of the 
BLR, making use of well sampled continuous light curves of 
the continuum and broad line emission, which are then cross- 
correlated. In Fig. [8] we compare our finding with variability 
data from different broad line QSO samples. The blue filled 
do ts are QSO variability data from the sample listed in Tab.6 
bv iKollats chnv. Zetzl & Die trichl (l2006bl) the red filled dots are 


data fromjR 


eterson et al 


12004 ) and Kasniet_ai]_| 2000h as listed 


m Tab.8 by KoliatschnwZetH & PietriclT i 2006h ). The compar¬ 
ison shows that the relation also holds for reverberation response 
of the broad line QSOs. With the additional literature data Fig.[^ 
now comprises a total of 61 sources. In the work we present here, 
however, the line response that we monitored in our small sam¬ 
ple is dominated by forbidden narrow line emission. Only the 
H/J line allows a separation between the broad and narrow line 
contribution. However, the median and mean derived from the 
BLR component only is in good agreement with the overall re¬ 
sult. We base our investigation on typically 2 measurements over 
a baseline in time of the order of 10 years (see Tab. lAlb . Re¬ 
verberation response of the NLR based on measurements over 
5-8 years has also been rep orted for the radio galaxy 3C390.3 
I Clayej_&_^^misteteil ll987ll and for the SI galaxy NGC 5548 
I Peterson et al.l2013h . Within the SDSS survey one can also find 
example of multi-epoch observations that are consistent with this 
picture. In Fie. lAlOl we show three epoch SDSS spectra for the 
NLSl sources J014412 and J022205. The variability estimates 
places them well into Fig. [8] 

Variability on time scale of years already sets an upper limit 
to the size of the corresponding region. Since the variability 
is rather significant, the line flux contribution of that region is 
also high. While there is a density and temperature gradient to¬ 
wards the nuclear position jOgle et alJ200(ll : TPenston et alJl99ol : 
iMorse et alj 19951) one also finds that the surface brightness dis¬ 
tribution of the NLR flux is very much c entrally peaked. Th e 
size of the NRL is luminosity denendant llSchmitt et alJfioO^ . 
however, the inspection of high angular resolution HST scans 
across S1 and S2 nuclei shows that the unresolved central peak 
may contain 50% or more of the m or e extended nuc l ear fo rbid- 
den line flux dBennert et alj|2006al l5). l^terson et ^ ( l2013l) find 
for NGC 554 8 that 94% of the N LR emission arises from within 
18p c (54 ly). QgleetalJ ll2000h derive a density profile of r“^ 
and IWalsh et alj i2008h find that the [S II] line ratio indicates a 
radial stratification in gas density, with a sharp increase within 
the inner 10-20 pc, in the majority of the Type 1 (broad-lined) 
objects. 

Hence, the main response to the continuum variation must 
therefore originate in a rather compact region with a diameter of 
the order of 10 light years. This is in good agr eement with the 
finding for 3C390..3 (Iciavel & Wamsteteill987h and NGC 5548 
dPeterson et al.l 12013 ). Since this formalism fully explains the 
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AL 

cont. 

Figure 8 . The continuum vaiiability plotted against the line variabil¬ 
ity. The solid and dashed black lines represent the relation = 

3 

6 (ALcon?.)- (see text). Here AL indicates the obsei'ved luminosity vaiia- 
tion in line and continuum. The filled colored dots are QSO variability 
data and continuum) from i Ko llatschny Zetzl & Dietrichl i20Q6bli . 
IPeterson et alj i2QQ4) . and iKaspi et alj j2QQ(l) . The Black crosses repre¬ 
sent the median values and their uncertainties as we derived them for 
our sample (^0 and other lines plus continuum) and as they are shown 
in Figs. |S]&|2] 


reverberation of the line emission with respect to the continuum 
variation (with a surprisingly close proportionality; see Fig. 
this suggests that the region illuminated by the nucleus is al¬ 
most fully spatially overlapping with the responding line emit¬ 
ting region. Hence, at larger distance the NLRs then show a 
much weaker response to the continuum variability most likely 
due to a combination of a lower overall luminosity and a lower 
volume filling factor. The close correlation to the non-stellar nu¬ 
clear continuum variability also implies that for the sources we 
investigated there are no other significant sources (i.e. star for¬ 
mation, extra nuclear shocks) but the nuclear radiation field that 
contributes substantially to the line emission of the reverberating 
region. 


7 CONCLUSIONS 

We investigated a sample of 18 sources. For 8 objects in 
Tab. [T] we performed a detailed analysis spectroscopy over the 
available optical wavelength range. For 10 sources we obtained 
suitable multi epoch line and continuum data from the litera¬ 
ture (see Tab. [8]and Figs.and O covering time scales from 5 
to 10 years. Fig. [8] demonstrates that our findings describe the 
variability characteristics of a total of 61 sources. From the line 
and continuum variability of these active galactic nuclei we find 
the following consistent picture that explains the differences be¬ 
tween the Seyfert 2/NLS 1 and the BLSl and QSO sample: The 
line luminosity L/,„e can be described as a reverberation response 
to the continuum luminosity Leant- In good agreement results of 
previous studies (references in section [631 > we find that the bulk 
of the NLR emission arises from within 10-20 light years and 
shows a reverberation response to variability of the nucleus. The 
differences in variability do not require a variation of the accre¬ 
tion rate between the two samples. It is mainly the differ¬ 
ence in black hole mass AM that is responsible for the difference 
between the samples. The increased variability of the line emis¬ 
sion can fully be explained by the dependency on the accretion 
rate ^M. 

at 
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Table 9. The line fitting results for J0938, J1203, 1158, J0911, J0802, J0354, J0153, and J0347. 


Sources Em. lines 

Ob. Wa. 

Inst. 1 


Inst.2 




Flux 

FWHM 

Flux 

FWHM 


A 

erg 

1 

A 

[10“^^ erg 

_ 1 2i 

A 





SDSS" 

LET 


J0938 

[O II] 23727 

4102 ± 1 

50.12 ±2.39 

5.66 ± 0.49 

43.08 ± 2.38 

6.05 ± 0.47 


[Ne III] 23868 

4258 ± 1 

11.28 ±0.91 

10.67 ± 0.93 

7.97 ± 1.14 

8.38 ±0.79 


Hy 2 4340 

4777 ± 1 

25.18 ± 1.58 

6.31 ±0.88 

18.95 ± 1.10 

6.91 ±0.71 


He II 24686 

5157 ±2 

8.50 ± 0.73 

16.81 ± 1.12 

4.04 ± 0.36 

10.66 ± 0.87 


HyS* 24861 

5350 ± 1 

50.28 ± 2.54 

5.61 ±0.49 

41.95 ±2.83 

5.91 ±0.73 


[O III] 25007 

5510 ± 1 

45.06 ± 1.57 

7.34 ± 0.57 

37.43 ± 2.01 

7.27 ± 0.66 


[Fe vii] 25721 

6298 ± 1 

4.60 ± 0.39 

10.56 ± 0.35 

6.55 ± 1.11 

11.74 ±0.64 


[Fe vii] 26087 

6699 ± 1 

11.55 ±0.75 

18.63 ±0.51 

11.19± 1.48 

19.65 ± 0.72 


[O i] 26300 

6935 ± 1 

4.94 ± 0.40 

6.75 ± 0.41 

5.30± 1.14 

9.14 ±0.83 


Hff 26563 

7224 ± 1 

208.83 ± 5.89 

7.75 ± 0.57 

164.56 ± 5.03 

6.43 ±0.61 


[N II] 26583 

7246 ± 1 

80.04 ± 2.57 

7.77 ± 0.52 

57.33 ± 2.70 

6.17 ±0.65 


[S II] 26716 

7393 ± 1 

28.90 ± 1.25 

7.46 ± 0.54 

23.08 ± 1.45 

5.92 ±0.52 


[Sii] 26731 

7409 ± 1 

19.28 ±0.96 

7.15 ±0.56 

16.48 ± 1.32 

6.02 ±0.51 




SDSS 

LET 


J1203 

[Ne v] 23345 

3898 ± 1 

20.42 ± 1.41 

5.75 ± 0.22 

27.82 ± 1.02 

5.18±0.31 


[Ne VI] 23425 

3993 ± 1 

55.54 ±2.84 

4.67 ± 0.25 

67.57 ± 1.93 

5.08 ± 0.24 


[On] 23727 

4344 ± 1 

83.03 ± 1.53 

6.14 ±0.14 

93.17 ±2.87 

6.19 ±0.22 


He 23969 

4625 ± 2 

38.36 ± 1.21 

4.92 ± 0.36 

48.01 ± 1.75 

4.98 ± 0.76 


Hy 24340 

5059 ± 1 

89.24 ± 1.49 

5.64 ±0.41 

81.84 ±2.04 

5.41 ±0.75 


HyS* 24861 

5665 ± 2 

154.54 ± 1.78 

5.88 ±0.75 

126.39 ± 1.11 

5.64 ± 0.97 


[0 III] 25007 

5835 ±2 

803.60 ± 6.39 

5.67 ±0.88 

681.90 ±7.66 

5.29 ± 0.85 


[0 i] 26300 

7343 ± 3 

12.53 ±0.75 

7.33 ±0.13 

11.25 ±0.54 

7.94 ±0.12 


Hof 26543 

7648 ± 2 

364.01 ± 3.60 

7.64 ± 0.55 

168.34 ±3.95 

7.12 ±0.71 


[N II] 26583 

7671 ± 1 

25.11 ± 1.63 

6.59 ±0.18 

13.17 ±0.94 

6.31 ±0.34 


[S II] 26716 

7827 ± 3 

17.09 ± 1.71 

7.00 ±0.11 

10.37 ±0.91 

7.40 ±0.14 


[S II] 26731 

7844 ± 3 

16.92 ± 1.28 

7.81 ±0.12 

8.98 ±0.88 

6.61 ±0.11 




SDSS 

LET 


J0911 

[Mg II] 22799 

4030 ± 1 

60.43 ± 1.53 

29.01 ± 0.95 

35.40 ±3.01 

29.01 ± 1.32 


[O II] 23727 

5372 ± 2 

11.80 ±0.84 

7.35 ± 0.48 

15.26 ±0.55 

7.77 ± 0.32 


Hy 24340 (broad) 

6258 ± 2 

15.71 ± 1.14 

54.14 ± 1.87 

18.42 ±0.77 

51.78 ±2.01 


Hy 24340 (nan'ow) 

6258 ± 2 

6.92 ± 0.75 

10.35 ± 0.69 

8.03 ± 0.41 

11.06 ±0.73 


H/? 24861 (broad) 

7007 ± 2 

46.05 ± 2.94 

54.14 ± 1.87 

52.55 ± 2.49 

51.78 ±2.01 


1^/14861 (narrow) 

7007 ± 2 

17.09 ± 1.06 

9.65 ± 0.62 

19.24 ±0.82 

10.59 ± 0.99 


[0 III] 25007 

7218 ± 1 

28.76 ± 1.53 

9.88 ± 0.84 

34.49 ± 1.55 

10.12 ± 1.02 


[0 i] 26300 

9082 ± 2 

4.75 ± 0.35 

20.22 ± 1.22 

5.43 ± 0.87 

16.14 ± 1.76 




SDSS 

LET 


J1158 

[O II] 23727 

5366 ± 1 

10.61 ±0.67 

8.35 ±0.51 

7.12 ±0.71 

9.85 ±0.88 


Hy 24340 (broad) 

6245 ± 3 

58.82 ±2.83 

75.32 ± 2.92 

29.63 ±2.16 

94.16 ±2.56 


Hy /14340 (narrow) 

6245 ± 3 

12.61 ±0.56 

14.83 ±0.81 

7.40 ±0.61 

14.35 ± 1.47 


H/? 24861 (broad) 

6994 ± 3 

201.01 ±3.31 

75.32 ± 2.92 

105.69 ± 4.46 

94.16 ±2.56 


/i4861(nan'ow) 

6994 ± 3 

28.21 ±0.96 

12.91 ±0.88 

17.05 ± 1.14 

13.98 ±0.94 


[0 III] 25007 

7208 ± 3 

59.41 ± 1.24 

12.01 ±0.91 

41.92 ± 1.88 

13.65 ± 1.18 


[O i] 26300 

9071 ± 3 

5.81 ±0.90 

6.55 ± 0.75 

3.98 ± 0.60 

8.49 ± 0.98 




SDSS 

LET 


J0802 

[O II] 23727 

6203 ± 1 

18.43 ± 1.23 

15.48 ± 1.31 

25.33 ±0.99 

12.91 ±0.89 


Hy 2 4340 (broad) 

7227 ± 2 

33.82 ±2.43 

80.11 ±2.02 

60.36 ± 0.95 

82.39 ± 2.89 


Hy 2 4340 (narrow) 

7229 ± 2 

8.62 ± 1.29 

12.00 ±0.79 

15.41 ±0.55 

11.29 ±0.47 


H/3 24861 (broad) 

8092 ± 1 

84.23 ±2.18 

80.11 ±2.02 

95.07 ± 2.76 

82.39 ± 2.89 


Hj6 24861 (narrow) 

8092 ± 1 

20.97 ± 1.63 

11.98 ±0.85 

43.22 ± 1.54 

13.65 ± 1.98 


[O III] 25007 

8331 ± 1 

16.40 ± 1.42 

11.51 ±0.81 

28.41 ± 1.62 

13.41 ± 1.12 




Beijing'’(2) 

LET 


J0153 

Hj3 24861 (broad) 

6453 ± 3 

420.80 ± — 

77.01 ± — 

299.92 ± 10.98 

51.36 ±3.98 


H^/14861 (narrow) 

6453 ± 3 

99.60 ± — 

18.34 ± — 

68.20 ±3.01 

14.96 ± 1.34 


[0 III] 25007 

6654 ± 1 

155.90 ± — 

13.37 ± — 

113.36 ±4.51 

14.12 ±2.45 


Continued on next page 


© 0000 RAS, MNRAS OQO.fTIfTTl 











Variability in Active Galaxies 


21 


Continued: The line fitting results for J0938, J1203, 1158, J0911, J0802, J0354, J0153, and J0347. 




Hiltner^(l) 

LET 


J0354 

[Ne v] T3346 

3544 ± 2 

26.00 ± 10.00 

4.51 ± 1.75 

11.91 ±0.97 

3.09 ± 0.68 


[0 III] 45007 

5183 ± 1 

413.00 ± 15.00 

5.51 ±0.94 

284.84 ± 4.42 

4.53 ± 0.26 


[O i] 46300 

6522 ± 1 

15.15 ± 1.34 

5.93 ±0.92 

8.70 ± 2.00 

6.60 ± 0.67 


[S II] 46716 

6952 ± 1 

14.45 ± 1.16 

5.82 ±0.34 

10.90 ± 1.00 

5.79 ± 0.49 


[S II] 46731 

6967 ± 1 

12.52 ± 0.96 

5.85 ±0.36 

8.40 ± 1.00 

5.09 ± 0.42 

J0347 



OHP'' 


LET 



Hy 4 4340 (broad) 

4475 ± 2 

752.63 ± 18.28 

51.78 ±2.68 

1094.56 ±9.19 

51.78 ±2.95 


Hy 4 4340 (narrow) 

4475 ± 2 

273.21 ± 5.23 

10.96 ± 1.23 

324.91 ±9.19 

10.01 ±0.55 


H/3 44861 (broad) 

5010 ± 1 

1644.70 ± 16.46 

51.78 ±3.54 

2280.33 ±21.87 

51.78 ±0.92 


/14861 (narrow) 

5010 ± 1 

1058.20 ± 11.89 

13.18 ± 1.02 

1398.38 ±20.46 

11.77 ±0.92 


[O III] 45007 

5156± 1 

951.10 ± 10.87 

12.71 ±2.33 

996.74 ± 20.80 

11.53 ±0.64 


[0 i] 46300 

6495 ± 2 

124.29 ± 2.67 

16.19± 1.11 

157.30 ±5.79 

18.78 ± 1.78 


References. (1 ) \GruDe. Pradhan &r Franii il200jl) . (2\ \Zhou et all 12002(1 . 
Notes: 


(*) Means HfS are narrow. 

(^) The Sloan Digital Sky Siirx’ey. 

(^) Beijing Obserx’atory. 

(^') Hiltner Telescope. 

Telescope of Obsen’atoire de Haute-Provence (OHP). 


Table Al. Broad and narrow emission lines variability. 


Type 

Sources 

Narrow line 
components 
line type Varib. % 

Broad line 
components 
line type Varib. % 

Narrow lines objects 

J0347 

Hy A 4340 

15.9 

Hy 2 4340 

31.2 



H/3 24861 

24.3 

1^24861 

27.8 


J1203 

Hy A 4340 

12.4 

Hy 2 4340 

_ 



Hy3 24861 

18.2 

24861 

- 


J0938 

Hy 2 4340 

24.7 

Hy 2 4340 

_ 



Hj3 24861 

16.5 

24861 

- 

Broad lines objects 

J1158 

Hy 2 4340 

31.8 

Hy 2 4340 

32.2 



Hj3 24861 

33.6 

24861 

31.4 


J0911 

Hy 2 4340 

19.8 

Hy 2 4340 

14.7 



HyS 24861 

11.1 

24861 

12.3 


J0802 

Hy 2 4340 

43.9 

Hy 2 4340 

60.5 



Hy3 24861 

51.4 

24861 

11.4 


J0153 

Hj3 24861 

31.5 

24861 

28.7 
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Figure Al. The same as Fig.[2but for the other galaxies. 
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Figure A2. Left: For a 0.8 arcsec slit we show for all sources the distribution of the intensity drops of the fainter exposures with respect to the brightest 
once. The data is shown in bin widths of 0.05 (fat dots) and for comparison half the value obtained with a bin width of 0.10 (black line). Right: 
Intensity drop due to the combination of variable seeing and slit offset. Case a indicates the level (dotted line) to which the spectra can be calibrated 
if for individual sources the fainter exposures are corrected to the level of the brightest once. For case /3 we assume that the brightest exposures all are 
subjected to the mean drop (fat black line) indicated by the statistics of the faintest with respect to the brightest exposures. The arrow and the dashed 
lines indicate the standai'd deviation from the mean. 
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Variability in Active Galaxies 
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Figure A3. The optical spectnam of J1203 and the result of calculating the difference and ratio of spectra from different epochs. 
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Figure A4. The optical spectrum of J1158 and the result of calculating the difference and ratio of spectra from different epochs. 
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Figure AS. The optical spectrum of J0911 and the result of calculating the difference and ratio of spectra from different epochs. 
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Figure A6. The optical spectrum of J0802 and the result of calculating the difference and ratio of spectra from different epochs. 
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Figure A7. The optical spectnam of J0354 and the result of calculating the difference and ratio of spectra from different epochs. 
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Figure A8. The optical spectnam of J0153 and the result of calculating the difference and ratio of spectra from different epochs. 
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Figure A9. The optical spectnam of J0347 and the result of calculating the difference and ratio of spectra from different epochs. 
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Figure AlO. Three epoch SDSS spectra for J014412 and J022205. These spectra allow us to derive vaiiability estimates of AcoHr.=0.18 and A/,=0.45 
for J014412 and Ac.o«r.=0.27 and A//„^,=0.26 for J022205. These values lie well above the low calibration uncertainties for the SDSS data (section [33^ . 
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